Mon. Not. R. Astron. Soc. 000, [HS] (0000) 



Printed 29 June 2010 



(MN style file v2.2) 



Radio emission from dark matter annihilation in the Large 
Magellanic Cloud 



o 

00 

6 



> 

in 

(N 
m 

\6 
o 
o 



Beatriz B. Siffert^*, Angelo Limone^, Enrico Borriello^'^, Giuseppe Longo^'^'^ 
and Gennaro Miele^'^ 

^ Dipartimento di Scienze Fisiche, Universita di Napoli Federico II, Via Cinthia, 9, 1-80126, Naples, Italy 
^INFN - Sezione di Napoli, Via Cinthia, 9, 1-80126, Naples, Italy 
^INAF-OACN, Via Moiariello, 16, 1-80128, Naples, Italy 

'^Department of Astronomy, California Institute of Technology, 1200 East California Blvd, Pasadena, CA 91125, USA 



29 June 2010 



ABSTRACT 

The Large Magellanic Cloud, at only 50 kpc away from us and known to be dark 
matter dominated, is clearly an interesting place where to search for dark matter 
annihilation signals. In this paper, we estimate the synchrotron emission due to WIMP 
annihilation in the halo of the LMC at two radio frequencies, 1.4 and 4.8 GHz, and 
compare it to the observed emission, in order to impose constraints in the WIMP 
mass vs. annihilation cross section plane. We use available Faraday rotation data from 
background sources to estimate the magnitude of the magnetic field in different regions 
of the LMC's disc, where we calculate the radio signal due to dark matter annihilation. 
We account for the e^e~ energy losses due to synchrotron. Inverse Compton Scattering 
and bremsstrahlung, using the observed hydrogen and dust temperature distribution 
on the LMC to estimate their efficiency. The extensive use of observations, allied with 
conservative choices adopted in all the steps of the calculation, allow us to obtain very 
realistic constraints. 
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1 INTRODUCTION 

Evidence for the existence of dark matter (DM) has been ob- 
served in various astrophysical systems, ranging from satel- 
lite dwarf galaxies to massive galactic clusters and cosmol- 
ogy. Possible cold dark matter candidates arise from theoret- 
ical models conceived to extend the Standard Model of ele- 
mentary particles and Weakly Interacting Mas sive P a rticles 
(WIMP) are the current main paradigm (see iFenel (|2010l ) 
for a recent review on dark matter candidates). 

The fact that direct observation of dark matter parti- 
cles has not yet been possible justifies all the efforts devoted 
to their indirect detection, i.e., the observation of anoma- 
lous components in the cosmic rays' spectru m that can be 
attrib uted to dark matter annihilation. See ICrocker et al 
2010t) 



Abdo et all (120101): ISergst rom. Fairbairn fc Pieril 
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Barger et al. (2002); Colafr ancesco. Profumo fc UUiol (|2006l . 



20071 ) for examples of such searches in the 7-ray, X-ray, radio 



and neutrino spectra. 

In fact, dark matter self-annihilation is expected to pro- 
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duce several Standard Model particles, among which elec- 
trons and positrons that, by interacting with the galactic 
magnetic field, emit synchrotron radiation in the radio band. 
We can estimate the intensity of the emission coming from 
a given direction in the sky by adopting models or using 
measurements to describe the following quantities: the dark 
matter density profile, which can be modelled by using the 
results of halo formation numerical simulations; the mag- 
netic field, estimated through techniques such as the anal- 
ysis of polarisation of radio and optical emission and rota- 
tion measures; and the interstellar radiation field and hydro- 
gen distribution, which are needed to account for the elec- 
tron/positron energy losses. The comparison of this emission 
with the observed radio emission then allows to impose con- 
straints on the values of the dark matter's mass, m-^, and 
its thermally averaged annihilation cross section, {gav). 

In order for this method to provide reliable constraints, 
the values of all the involved parameters need to be known 
very well, which often does not happen. Furthermore, ide- 
ally, the comparison between the theoretical result and the 
observed emission should be made after all the known as- 
trophysical foreground has been subtracted from the obser- 
vations. In practice, however, such subtraction is subject to 



2 Beatriz B. Siffert et al. 



several uncertainties, since often we cannot accurately model 
all the components contributing to the foreground. 

In addition, the method presents an intrinsic shortcom- 
ing: we can only calculate the contributions from DM an- 
nihilation integrated along the line of sight, so all the DM 
signal may be 'averaged' away by low emitting regions that 
happen to fall along the same direction in the sky. 

In the present paper, we present constraints obtained by 
applying this method to the Large Magellanic Cloud (LMC) 
and comparing the results with radio observations at 1.4 and 
4.8 GHz. Unless not otherwise stated, we will follow the for- 
malism described in iBorriello. Cuoco &: Mield { 20091 ) in our 
calculations. We analyse two possible WIMP annihilation 
channels, a hadronic and a leptonic one, the latter having 
been recently proposed as a possible cause of the anomalous 
cosmic ray electron / positron spectrum measured by Pam ela 
l|Adriani et al.ll2009l) and Fermi-LAT (|Abdo et al.il2009l '). 

By extensive use of available observations of the LMC 
in different frequencies, we are able to obtain realistic infor- 
mation to describe all the inputs necessary to estimate the 
DM annihilation signal. In this way, we escape from the 'too 
many hypotheses' problem that has recurrently appeared in 
the existing literature on the subject. When we do have to 
make a hypothesis, we follow the most conservative path, 
making sure that our choice will not overestimate the sig- 
nal. 

The LM C, at a distance of only ~ 50 kpc from us 
(|Alvesll2004l '). is probably one of the best studied galax- 
ies in almost all frequency bands. Rotation curve data 
imply that the LMC total mass inside a 9 kpc radius 



is ~ 1.3 X 
neutral gas 



10'" Mq, 
masses is 



while the sum of the stellar 
estima ted to be ~ 3.2 x 10^ 
Beslal bOOgl ). Therefore, 



and 
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the 

disc 
sky 



(|van der Marel. Kallivavalil 
LMC must be dark matter dominated. 

The fact that the LMC is nearly face-on (its 
forms an angle of ~ 35° with the plane of the 
llvan der Marel fc Cionill200ll )). minimises the line-of-sight 
integration problem mentioned above, since, in this case, 
the integration spans only over the thickness of the galaxy. 
A nearly face-on view also allows us to clearly distinguish 
between high emitting and low emitting regions within the 
disc, i.e. between regions where different astrophysical pro- 
cesses are taking place and different astrophysical fore- 

ground s exist. 

In IBorriello etahl (|2010l '). some of the authors of the 
present paper looked for DM annihilation from low radio 
emitting regions (called 'radio cavities') within the nearly 
face-on Local Group member Messier 33 (inclination 56°). 
Following this reference, a good candidate radio cavity 
should present low radio emission relative to the rest of the 
disc and should not be very far from the centre of the galaxy. 
While the first criterion guarantees that the astrophysical 
foreground in that region is relatively small, the second one 
is required so that the DM density is non negligible and that 
the possibility of a high magnetic field is not ruled out. Al- 
though this study produced constraints com parable to those 
obtain ed for the Milky Way (see for example IBorriello et al.l 
Hooi)), the lack of knowledge of the values of the astrophys- 
ical environment governing the diffusion of the electrons and 
positrons, in particular the magnetic field, introduces several 
uncertainties on the results. 

In the current analysis, in contrast, we use recent accu- 



rate rotation mea sures (RMs) from poin t sources behind the 
LMC obtained in lGaensler et af] (|2005l ), through which we 
can directly estimate the magnitude of the magnetic field 
in several regions for the disc of the LMC. We focus our 
analysis on 23 regions on the LMC's disc with galactocen- 
tric distances < 8 kpc for which a KM has been measured. 
Although these regions are not chosen to be radio cavities, 
the fact that a background radio source has been seen there, 
guarantees that they present low radio emission in compar- 
ison to other parts of the galaxy (intensity at 1.4 and 4.8 
GHz < 2 mjy/beam). For each of these regions we calculate 
the expected DM emission at 1.4 and 4.8 GHz and com- 
pare the results with the observed values without performing 
any foreground subtraction (which is a conservative choice, 
since removing part of the observed emission would always 
strengthen the constraints). 

Previous searches for dark matter annihilation 
signa ls from the LMC have analy s ed b oth the 7- 
rav JCondold 1 19941: iTasitsiomi et all 120041 ) and radio 
(|Tasitsiomi. Gaskins fc Olintd 2'o04h bands. In particular, 
in the latter the flux density in several radio frequencies due 
to annihilation into the hadronic channel integrated over 
the enti re galaxy is calcu l ated. However, unlike what is done 
here, in ITasitsiomi et al] (|2004l ) a simple constant magnetic 
field inside the entire volume of the galaxy is assumed. 
Furthermore, electron/positron energy losses due to Inverse 
Compton Scatte r ing a nd bremsstrahlung are neglected in 
ITasitsiomi et al.l (|2004l ) although, as we show in section |4j 
these processes can be very important in several regions of 
the LMC. When appropriate, we will comment on the effect 
of these assumptions and compare our results with theirs. 

The paper is organised as follows. In section [2] we 
present the radio data used in our analysis and the 23 re- 
gions studied. Section O is dedicated to the distribution of 
dark matter in the halo of the LMC and presents the results 
of fits to different dark matter profiles. In section U we dis- 
cuss the energy loss processes important to our calculations. 
In section [5] we present our results and discuss them and 
compare with previous ones in section [B] 



2 RADIO DATA AND SELECTION OF 
REGIONS 

We used radio observations of the LMC at two different fre- 
quencies, 1.4 and 4.8 GHz. At 1.4 GHz, we used the mosaic 
image obtained with the Australia Tel escope C ompact Ar- 
ray (ATCA) and the Parkes Telescope (iHughes e t al. 2007i), 
shown in the left panel of Fig. [1] The HPBW for this image 
is 40 arcsec, which corresponds to ~ 10 pc at the LMC's 
position, and its sensitivity is estimated to be between 0.25 
and 0.35 mJy/beam. At 4.8 GHz, we used the observations 
made with the ATCA described in Dickel et_al. (2005) and 
shown in the right panel of Fig. [T] The HPBW is 33 arcsec 
and the estimated sensitivity is 0.28 mJy/beam. 

We estimated the radio emission due to dark matter an- 
nihilation at these two frequencies coming from 23 circular 
regions with radius 180 arcsec within the LMC's disc. The 
selection of regions was made according to the results pre- 
sented in lGaensler et all l|2005l ). where the Faraday rotation 
for 291 polarised sources behind the LMC was measured, 
100 of which happen to lay directly behind the LMC. Af- 
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Figure 1. Radio maps in logaritlimic scale of tiie LMC at 1.4 GHz l|Hughes et al.1120071 ') and 4.8 GHz jPiclcel et al.ll2005h . Tlie numbers 
indicate the positions of the 23 selected regions where the Faraday rotation of background sources is measured. The label 'GC' indicates 
the assumed LMC's kinematic centre. 



ter subtracting the effects due to the foreground Faraday 
rotation in the Mi lky Way, a map of ro tation measures is 
obtained (fig. 1 in iGaensler et alT ()2005l ll. where circles of 
various sizes represent the magnitude of the RM at each 
position. We chose the positions of 23 circles to apply our 
method, selecting them according to their diameter, which 
is proportional to the RM at that position (large RM im- 
plies large magnetic field and therefore large synchrotron 
emission), and distance from the centre of the LMC (as we 
move away from the centre the density of DM decreases). 
In our calculations we used the LMC's kinematic centre, 
with coordinates RA= 05 h 17 m 36 s and DEC= -69°02' 
l|Kim et al.111993 ). 

Since the data published in ICaensler et al.l (|2005l ) are 
not yet of public domain, we had to obtain the positions and 
the RM of each of the 23 regions directly from the map in 
their fig. 1, using the fact that the maximum RM measured 
is +247 ± 13 rad m~^. To help us determining the coordi- 
n ates of each region, we used th e list of sources published 
in lMarx. Dickey fc Meboldl ([l993), which contains 113 com- 
pact radio sources detected with the ATCA at 1.4 GHz in 
and behind the LMC, in a region comprehended between 
05hl7m < RA < 05h51m and -71° < DEC < -67°. It is 
estimated that among these 113 sources, 15 are in the LMC, 
so most of them are background objects. 

Once RMs have been estimated for each region, the cor- 
responding magnetic field was obtained using the definition: 

RM = k J Brcgii ds, (1) 

where k = 0.81 rad m~^pc~^cm^ fj,G~^, Ue is the electron 



number density, Brcg|| is the regular component of the mag- 
netic field parallel to line-of-sight and s is the distance along 
the line-of-sight. 

Assuming that the regular component of the magnetic 
field, Brcg, hes solely on the plane of the disc, we conclude 
that Brcgii is produced solely by the galaxy's inclination. 
The mean value of Brog, weighted according to the electron 
density Ue along the line-of-sight is then given by 

RM 



-Bros — 



k T>M sini ' 



(2) 



where VJ^ is the dispersion measure, defined as J Ue ds and 
i ~ 35° is the LMC's disc inclination. In our calculations 
we used the mean total dispersion measure calculated from 
five ra dio pulsars in the LMC obtained in ICrawford et al.l 
ll200ll), "DM = 100 cm ^ pc, following what was done in 

iGaensler et al] (|2005|) 

Inl Gaensler et al.l (|2005l ) . the random component of the 
LMC's magnetic field, Bran, was estimated to be ~ 3.6 Brog, 
so that Btot = \/ Bre g + Bran — 3.7Brcg. As mentioned in 
IGaensler et"all (120051 ) . in regions where the magnetic field 
and rie are correlated, the above calculation overestimates 
Brog by a factor 2 and underestimates Bran by the same 
factor. 

The 23 selected regions are shown in Fig. [1] super- 
imposed to the radio maps. They are numbered from the 
smallest to the largest galactocentric distance. In Table[T]we 
list for each region: its equatorial coordinates; correspond- 
ing ba ckground source name as in [m arx. Dickev fc Meboldl 
(|l997l l. if any; distance to the centre of the LMC along the 
plane of the disc, assuming a distance to the LMC of 50 
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Table 1. The 23 regions studied. Intensity is measured in mjy/beam. 





RA (J2000) 


DEC {J2000) 


Name 


d (kpc) 


Brcg (mG) 


Btot (/iG) 


U^id'^leV cm-») 


nn (cm ^) 


Il.4 GHz 


I4.8 GHz 


1 


05:17:43.541 


-69:35:17.07 


MDM3 


0.56 


2.45 


9.14 


1.508 


0.07 


1.83 


1.58 


2 


05:13:14.664 


-69:31:17.19 




1.21 


2.76 


10.33 


1.346 


0.14 


1.67 


1.27 


3 


05:17:22.087 


-70:20:41.31 


MDMl 


1.33 


1.59 


5.96 


1.098 


0.07 


1.56 


1.01 


4 


05:18:51.899 


-67:45:54.21 


MDM8 


1.36 


3.61 


13.51 


1.030 


0.07 


0.94 


0.67 


5 


05:23:46.213 


-68:45:33.88 


MDM20 


1.51 


1.70 


6.36 


1.207 


0.14 


2.18 


1.47 


6 


05:10:45.796 


-68:05:11.09 




1.71 


2.23 


8.34 


1.135 


0.07 


1.40 


0.69 


7 


05:21:50.892 


-70:35:46.96 


MDM18 


1.71 


2.45 


9.14 


0.439 


0.07 


0.84 


0.43 


8 


05:11:01.258 


-69:34:40.98 




1.72 


2.66 


9.93 


0.857 


0.14 


1.06 


0.60 


9 


05:22:57.600 


-67:29:16.80 


MDM17 


2.17 


5.32 


19.86 


0.921 


0.07 


0.80 


0.48 


10 


05:10:21.596 


-67:17:12.72 




2.21 


3.08 


11.52 


0.997 


0.54 


1.52 


0.58 


11 


05:27:09.571 


-67:49:06.69 


MDM30 


2.73 


2.34 


8.74 


1.266 


1.93 


2.09 


1.69 


12 


05:28:40.638 


-70:35:34.72 


MDM32 


2.78 


2.23 


8.34 


1.233 


0.14 


1.50 


1.23 


13 


05:06:19.726 


-71:08:26.64 




3.67 


4.04 


15.10 


0.637 


0.07 


0.63 


0.28* 


14 


05:26:09.659 


-66:00:02.93 




3.95 


4.25 


15.89 


0.569 


0.42 


1.63 


1.73 


15 


05:02:44.925 


-66:00:28.26 




4.19 


4.57 


17.08 


0.643 





0.68 


0.36 


16 


05:34:13.644 


-67:55:04.42 


MDM58 


4.23 


2.23 


8.34 


1.106 


0.14 


2.11 


1.38 


17 


04:56:57.832 


-68:50:35.85 




4.79 


1.70 


6.36 


0.985 


0.24 


1.00 


0.78 


18 


05:41:11.743 


-68:03:33.24 


MDM81 


5.77 


3.72 


13.90 


0.601 


0.14 


0.78 


0.38 


19 


05:40:54.847 


-67:41:04.43 


MDM82 


5.85 


3.72 


13.90 


0.694 


0.14 


1.12 


0.53 


20 


04:51:36.413 


-66:51:26.34 




6.07 


2.34 


8.74 


0.900 





0.64 


0.28* 


21 


05:45:23.757 


-68:42:57.08 


MDM95 


6.56 


2.76 


10.33 


0.878 


0.14 


1.35 


0.53 


22 


05:50:19.691 


-69:00:39.60 




7.64 


2.45 


9.14 


0.494 


0.14 


1.00 


0.62 


23 


05:51:58.795 


-69:45:06.22 


MDM113 


7.92 


2.45 


9.14 


0.461 


0.14 


0.46 


0.28* 



kpc; regular and total magnetic fields estimated from the 
RM; interstellar radiation field on the disc (see section |4]); 
the atomic hydrogen number density (see section and in- 
tensity measured at 1.4 and 4.8 GHz. The values of intensity 
at 4.8 GHz marked with an asterisk were lower than the cor- 
responding sensitivity and were set to the sensitivity value, 
0.28 mJy/beam. 

To model the behaviour of the magnetic field along the 
z-axis (perpendicular to the disc plane) we use an exponen- 
tial decay of the following form: 

B(z) = Bdi.ce-l^l/"°, (3) 

where zq is the scale height of the magnetic field, which 
can be up to 4 times larger than the scale height of the 
thick synchrotron disc (in the case of equipartition between 
cosmic ray s and magn etic fields and a synchrotron spectral 
index ~ 1) l|Beckll200j ). The thick disc has typically a scale 
height of ~ 1.8 kpc, evidenced from its direc t measurement 
in edge-on galaxies (|Dumke fc Krause]|l998l ). i n agreemen t 
with its value for the Milky Way, ~ 1.5 kpc l|Beckll200ll ). 
Therefore, the magnetic field scale height can be as large 
as 6-7 kpc. Since the LMC is a smaller galaxy (with radius 
~ 3 times smaller than the Milky Way's), we use in our 
calculations the value zo = 2 kpc. 

Let us notice that the exponential decay in eq. ([Sjl is 
a conservative choice and can only underestimate the value 
of the field, since the value of the magnetic field obtained 
through the RM is already mediated along the line-of-sight. 



3 DARK MATTER DENSITY PROFILE 

The intensity of the radio emission due to WIMPs annihi- 
lation coming from a given region within the LMC's disc 
depends on how these particles are distributed in the halo 
of the LMC. In fact, if p{r) is the density of WIMPs as a 
function of galactocentric distance r, the emissivity at a give 



frequency will be proportional to p^{r). We need, therefore, 
to estimate p(r) in order to estimate the emission due to 

DM an nihilation. 

In 'Tasitsi omi et al.l l|2004h . different DM density pro- 
files are studied using measurements of the LMC's rota- 
tional velocity field. Follo wing this appro ach, we used the 
result s obtained from HI (|Kim et al.lll99^ ) and carbon star 
data (|Alves fc NelsonI I2OO0I ) to estimate p(r). The HI ve- 
locity field is composed of 26 data points at galactocentric 
distances between ~ 0.05 kpc and ~ 3 kpc. The carbon stars 
obtained from 422 stars results in 4 data points, with galac- 
tocentric distances between 4.0 kpc and 8.2 kpc. We fit to 
this data set six different DM density profiles suggested in 
the literature: 

• A simple isothermal sphere cored profile: 

• the Navarro-Frenk-White (NFW) profile, derived 
from N-body simulations of CDM haloes for struc- 
tures ranging from dwar f gala x ies to clusters of galaxies 
(jNavarro. Frenk fc Whitd 1 19951 . Iigod l , which diverges as 
p oc in the central region of the galaxy: 



pNFw(r) = 




• the Moore et al. profile (|Moore et al.l Il999l ). derived 
from independent higher resolution CDM halo simulations, 
which has a steeper slope than the NFW profile for the cen- 
tral regions of the galaxy (p oc r~^'^): 



PM(r) = 




• the Burkert profile (|Burkertlll995f ). i.e. a cored profile 
derived from rotation curve fits of dwarf spiral galaxies: 
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Figure 2. Velocity field data points superimposed with the re- 
sulting fitted curves for the isothermal sphere, NFW, Burkert and 
Hayashi et al. density profiles. 
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Figure 3. Results of fits for the Einasto density profile for a ■ 
0.075, which resulted in the smallest value of x^/d-o.f., and a ■ 
0.3. 



PB(r) = 



Po 



• the Hayashi et al. profile (jHavashi et al.|[2003l '). which is 
a modified NFW profile derived from N-body simulations of 
the evolution of CDM subhaloes undergoing tidal stripping 
while orbiting around a larger halo: 

ft 



PH{r) 



PNFw(r) 



Po 



1 + 



rg > 



1 + 



where po = ft. Po.nfw; 

• and, finally, the Einasto profile, with three free para m- 
eters, first introduced in 1965 by Einasto (|Einastclll965l ') to 
describe the distribution of stars in our galaxy, and later 
used to describe the result s of N-body simulations of DM 
haloes (jNavarro et al.ll2004 ): 



PE(r) = poexp<^ - 
a 



- -1 



The NFW, the Moore et al. and the Hayashi et al. profiles 
were already tested in iTasitsiomi et ah 

The values of the free parameters were determined by 
fitting each profile to the rotation velocity field using that 
the velocity V{r) at an orbit of radius r, inside of which a 
DM mass M(r) is enclosed, is given by: 



V{r) = 



GM{r) 



(4) 



The Moore et al. profile gave the worst fit to the data, 
with a x^/d-o.f. ~ 8. The fitted curves for the other profiles 
can be seen in Figs. [2] and [3] and the resulting values for 
the free parameters for each density profile are displayed in 
Table [2l as well as the value of the corresponding /d.o.f. 
For the purpose of fitting the Hayashi et al. profile, we fixed 
ro,NFW ~ 9.04 kpc, its best-fitting value previously found 
for the NFW profile. 

We can see that both the isothermal sphere and the 
Burkert profile did not fit well the data, resulting in 
xVd.o.f. > 4. 



Table 2. Values of the parameters po and tq and jd.o.f. ob- 
tained from the fit to rotation curve data for each profile. 



Profile 


PO (xlO'' M0kpc-^) 


ro (kpc) 


X'id.o.f. 


NFW 


8.18 ± 2.67 


9.04 ± 2.43 


1.26 


Iso. Sphere 


326 ± 65 


0.52 ± 0.076 


4.23 


Burkert 


289 ± 51 


1.06 ± 0.13 


4.54 


Hayashi et al. 


8.16 ± 0.30 


6.36 ± 2.08 


1.17 


Einasto 


0.0015 ± 0.0008 


824 ± 381 


1.52 



For the Einasto profile, we found that good fits (with 
values oi I d.o.j . around 1.5) can be obtained for values of 
Q between 0.06 and 0.1. The smallest value of /d.o.f. was 
obtained for a — 0.075 and is displayed in Tabled We can 
see in Fig.O that the resulting curve for a = 0.075 fits very 
well the rotation curve data for small values of the radius, 
but does not reproduce its behaviour at larger galactocentric 
distances. We also plot in Fig. [3] the curve obtained with 
a — 0.3. In this case, the fit results in pE ~ 6.38 ± 1.32 x 10® 
M0 kpc"^ and te = 4.16 ± 0.67 kpc. Although this curve 
reproduces better the global behaviour of the data, we get 
a larger value of x^ /d.o.f. = 2.05. This happens because 
the small radius points (with smaller error bars) are better 
fitted by the curve with a = 0.075. 

In any case, the NFW and the Hayashi et al. profiles 
fit the rotation curve better than the Einasto profile. The 
behaviour of these two density profiles (with the values of 
the free parameters set to the ones in Table [21) as a function 
of the galactocentric distance r is shown in Fig. [4] We can 
see that the DM densities predicted by the two profiles are 
practically the same for r < 2 kpc. 

Nevertheless, at the largest distance sampled by our 
background sources, 7.92 kpc, the difference between the two 
density profiles is much smaller than one order of magnitude: 
PNFw(7.92 kpc) ~ 2.8 pH (7.92 kpc). 

We therefore chose to use the Hayashi et al. density 
profile in our calculations, both because its fit to the data 
resulted in the lowest value of x^ /d.o.f. and because it pro- 
vides a DM density value lower than the one provided by 
the NFW profile for most of the regions we are interested in 
(conservative choice). 



6 Beatriz B. Sijjert et al. 



1x10* 



7 2x10' 

O 



1 X 10' 
5x 10'' 

2x 10'' 
1x10'' 



NFW 



Hayashi et al. 



1 2 3 4 5 6 7 8 

r (kpc) 

Figure 4. NFW and Hayashi et al. density profiles as a function 
of r for the best-fitting parameters. The dots show the positions 
of all the 23 regions. 



4 ENERGY LOSS PROCESSES 

While diffusing away from the production site, several en- 
ergy loss processes act on the e^e~ fluid. Our calculations 
accounted only for the three fastest energy loss processes: 
the synchrotron emission, the inverse Compton scattering 
(ICS) off the background photons and bremsstrahlung. We 
do not include Coulomb energy losses. 



4.1 Synchrotron 

An electron or a positron with energy Ee diffusing in a region 
filled with magnetic field B will lose energy via synchrotron 
radiation with the following characteristic time-scale: 



B 



ICeVi ' 



where r°yn = 3.95 x 10^^ 



4.2 Inverse Compton Scattering 

Similarly, the characteristic time for the ICS is given by 

Ureui 



TICS = Tics 



eVcm 



[GeVj ' 



where ti'cs ~ 9.82 x 10^^ s and Umd is the energy density of 
the Interstellar Radiation Field (ISRF). 

In order to estimate the intensity of the ISRF at the 
position of each region, we used the map of dust t empe rature 
(Td) in the LMC presented in iBernard etaH (|2008l '). The 
energy density of the ISRF is proportional to T^^^ . We used 
P = 2 iBoulanger et al]| 19961 ) . To obtain the proportionality 
constant, we use the fact that the dust temperature in the 
solar neighbourhood is 1 7.5 K and the local value o f the 
ISRF is 0.539 eV cm~ ^ (|Weingartner fc Drainel 120011 ). We 



can therefore estimate the the ISRF at the position of each 
one of the 23 regions we are interested in. 

Following the same conservative approach we adopted 
in the case of the magnetic field, we impose that the value 
of the ISRF obtained in this way corresponds to its value 
on the disc, , and assume that, for each region, Umd 
decays exponentially with z: 



(5) 

where ho is the disc scale height, which is ac tually a function 
that grows with r and can be modelled as (jAlves fc NelsonI 
[2000): 



ho{r) = /io(0) e'/^ kpc. 



(6) 



where ho{0) — 0.14 kpc is the value of the disc scale height 
in the centre of the LMC and ^ = 2.24 kpc. 

Table [U shows the values of U^^S" for all the 23 regions 
under study. 

4.3 Bremsstrahlung 

The char acteristic time for the bremsstrahlung process is 
given by (|Longaiij|l997h : 



Tbr, 



'Vcm~3y ' 



where r^,.^^ = 1.17 x 10^^ s and tih is the hydrogen number 
density. 

To estimate the value of nu in our selected regions, 
we used the tabl e of HI clouds in the LMC presented in 
iKim et al.l (|2003), which contains the coordinates, radius 
and mass of each cloud. We assumed that the clouds are 
spherical and estimated their density by simply dividing 
their mass by the volume of a sphere with the reported ra- 
dius. We then checked if each one of the 23 regions fell inside 
a cloud and, if so, associated to that region the density of 
the cloud. If a region fell inside more than one cloud, we 
associated to that region the density of the highest density 
cloud, which is a conservative choice. Regions that did not 
fall inside any cloud were assigned with nu ~ 0. The values 
of TiH for each region obtained in this way are reported in 
Table [T] 

We neglect the contribution of ionised hydrogen in our 
calculations. This is a reasonable approxi mation in view 
of the depolarisation results obtained in ICaensler et al.l 
which indicate that there are no bright individual 
HII regions in the directions where the RMs were observed. 

To include the bremsstrahlung in the calculat ions, we 
need to modify the ^(r) function defined in lBorrieUo et al.l 
(|2009h to: 



M(r) 



syn t-^rad (r) ^ T^ 





syn 



T°cs eVcm ^ 
-1 



i/Q B{r) nu 
V fiG cm~' 



+ 



where v is the frequency and uq = 3.7 x 10^ Hz. 

Fig. [5] compares the energy loss rate, b = Ee/r, for 
each of the three energy loss processes, at the position of 
all the 23 regions. For each region, we plot the values of the 
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Figure 5. Relative contribution of each energy loss process calculated at the position of each region at 1.4 GHz (see text). The darkest 
and lightest colours (bottom and upper parts) show the synchrotron and the bremsstrahlung contributions, respectively. The middle 
parts show the ICS contribution. 



normalised energy loss rate for the synchrotron (calculated 
using -Breg), ICS and the bremsstrahlung processes. All the 
values were evaluated on the plane of the disc {z = 0) and we 
used the frequency peak approximation Ee = -y/ Brag x v/uo 
GeV, where Brog is measured in /iG. Since the results were 
very similar for v = 1.4 and v = 4.8 GHz, we show only the 
results obtained for 1.4 GHz. 

While ICS is clearly the dominant process in all regions 
we are studying, the contribution of bremsstrahlung seems 
to be the the smallest one. It is however, surely not negli- 
gible in at least 3 regions (10, 11 and 14). These processes 
therefore, cannot be neglected. 

If we remember that the regions are numbered accord- 
ing to their galactocentric distance, we immediately notice 
in Fig.[5]the apparent lack of correlation between the energy 
loss rates (and therefore of the magnetic field, the ISRF and 
the hydrogen distribution) and the galactocentric distance, 
which is clearly a consequence of the irregular nature of the 
LMC. 



5 RESULTS 

We considered two possible WIMPs annihilation channels: 
XX ~^ bb, in which electrons and positrons will be produced 
by decaying muons (/it~ — > e~ VeVfj.) and anti-muons (/i+ 
e^VuVe) produced in pions decays (tt" — > ''■nd vr"^ — > 

H'^Vfi) and the leptophilic channel XX ~^ ■ 

Leptophilic channels have recently raised interest in 
view of the experimental results on the electron/positron 
cosmic ray spectra. While Pamela observ ed an unexpected 
rise in the positron fraction ( |Adriani et al. 2009), Fermi- 
LAT observes a de viation from a simple power-law spectrum 
l|Abdo et al.ll2009|), thus confirming the previous results ob- 
tained by HESS (jAharonian et al1l2009l ). If these results are 
to be interpreted as due to DM annihilation in the galac- 
tic halo, one needs to consider leptophilic channels and high 
mass scales. 

The synchrotron intensity at a frequency u due to DM 
annihilation coming from a region inside a solid angle dfl on 
the LMC's disc is: 



dl cosi f . , , , , , 

where i is the disc inclination, s is the distance along the 
line-of-sight and jv {r, z) is the synchrotron emissivity at a 
position in the LMC's halo with galactocentric distance r 
along the disc and height z above or below the disc. The 
term cosi accounts for the fact that the line-of-sight is not 
parallel to z. 

The expre s sion o f the emissivity is derived in detail in 
iBorriello et all (|2009( ). 

We evaluated eq. (O for all 23 regions using v = lA and 
4.8 GHz. The most constraining results were obtained for 1.4 
GHz and can be seen in Fig. [6] for both annihilation chan- 
nels considered. Fo r comparison, we show also the best con- 
straint obtained in IBorriello et all l|2010l ) for M33 using an 
NFW profile and assuming equipartition between magnetic 
fields and cosmic ray s, and the best constraint obtained in 
IBorriello erahl (|2009l ) for the Milky Way using the xX ^ bh 
channel. For the xX ~^ M^A*" channel we show the best 
constraint obtained for the Milky W ay using the same for - 
malism and observations described in IBorriello et all (|2009l ) 
and also the favoured region obtained when one attributes 
to DM annihilation the experimental results described above 
(jMeade et al.ll201Gl ). 



6 SUMMARY AND DISCUSSION 

We have imposed constraints on the m^-{oAv) plane using 
radio observations at 1.4 GHz and 4.8 GHz of the LMC 
and analysing two different DM annihilation channels, a 
hadronic and a leptonic one. 

The existence of high resolution observations of the 
LMC in several frequency bands has allowed us to obtain 
most of the information needed to calculate the DM annihi- 
lation signal, making the least possible number of hypothe- 
ses in all the steps of the calculation. Being able to escape 
from this problem and, when necessary, making conserva- 
tive assumptions, we have, therefore, obtained very robust 
results. 

In all cases studied, the best constraints were obtained 
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Figure 6. Constraints obtained with tiic XX ^ channel (left) and the XX ~^ M^M" channel (right). The solid lines indicate the 
best constraints obtained with the LMC using either Stot or -Brcg- The dashed lines indicat e the best con straints obtained with M33 
llBorrie llo ct al. 2010) and the dot-dashed ones indicate those obtained with the Milky Way |Borriello et al. 2009). The s haded regions 
are t he forbidden ones. The ellipse shaped region indicates the regions favoured by the Pamela/HESS/Fermi-LAT results l lMeade et al.l 
l20ld) . 
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Figure 7. Same as Fig. |6] using the isothermal sphere DM density profile with the parameters obtained in lAlves fc Neisoiil l|200Ch . 



using 1.4 GHz, since the number of e~ and e produced by 
DM annihilation decreases with energy. For higher frequen- 
cies to produce competing constraints, the observed inten- 
sity at those frequencies should be much lower than at 1.4 
GHz. 

From Fig. [6l we can see that the constraints on the 
m^-{aAv) plane imposed from the analysis of the LMC are 



stronger than the ones obtained with M33 and with the 
Milky Way over the most of mass range considered when, 
excluding very low mass regions. 

In iTasitsiomi et al.l (|2004l ) the estimated DM annihi- 
lation signal at different frequencies was obtained, fixing 

= 50 GeV and {aAv) = 2 x 10"^® cm^s~^ and using 
a constant magnetic field of intensity either B = 5 nG or 
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B = 18.4 /iG. We can estimate from their fig. 5 that at 1.4 
GHz they obtain (aAv) < 4 x 10"^^ cm^s"^ for B = 5/iG 
and (aAv) < 3 x lO"^'^ cm^s"^ for B = 18.4 /iG. From 
our Fig. [6] we can see that for = 50 GeV our best con- 
straints using the hadronic channel are (aAv) < 9.25 x 10 
cm^s~^ using only B^eg and {aAv) < 1.10 x 10~^* cm^ s 
using Btot . Th is result can be explai ned by the hypothe- 
ses adopted in ^Tasitsi omi et al.l (|2004l ) , which have overes- 
timated the DM annihilation signal. In the first place, in- 
spection of Table [1] shows that imposing a constant mag- 
netic field of 18.4 fj. G for the entire volume of the LMC is 
not realistic. Even if we take into account only the Btot 
column, we see that more than half of the regions present 
Btot < 10 /i G. As for the hypothesis S = 5 /i G, we can see 
from Fig. [5] that the bremsstrahlung and, in particular, the 
ICS processes cannot be neglected (see for example regions 
9 and 15, for which Brog = 5.32 /i G and 4.57 fi G, respec- 
tively). In fact, if we do neglect these two processes in our 
calculations, our results improve by a factor ~ 3. 

Finally, we would like to comment on one possible short- 
coming of our analysis. When fitting the LMC's rotation 
curve, we took into consideration only the DM mass, ignor- 
ing the luminous matter. To understand the effect of this 
ch oice, we redid our calc ulations using the results obtained 
in lAlves fc NelsonI (|2000l ). where an isothermal sphere den- 
sity profile is fitted to the LMC's rotation curve after sub- 
tracting the disc component (stars + gas). The resulting 
values of the free parameters are piso ~ lO^M© kpc~'^ and 
riso ~ 1 kpc. Using this density profile, the best constraints 
were obtained at 1.4 GHz and can be seen in Figs. 7(a 



and 7(b) 



We can see that, accounting for the luminous matter 
in this way, our results are not very much affected and the 
conclusions previously described are still valid. 
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